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SUMMARY 


"■'^This  memorandum  describes  the  current  version  of  the 
program  used  for  6  component  measures snts  using  strain  gauge  balances 
and  the  PDP8/I  installation  in  the  Transonic  Wind  Tunnel.  The 
program  is  an  updated  version  of  ARL  Tech.  Memo.  A264^and  has 
been  in  use  for  several  years.  Changes  to  the  original  program 
include  on-line  display,  model  weight  correction,  and  operation 
with  the  Pollock  strain  gauge  equipment*^-  Although  the  8/1  is 
scheduled  for  replacement,  it  may  be  some  time  before  the  change¬ 
over  to  a  new  installation  can  occur.  Also,  some  of  the  program 
detail  may  be  useful  to  the  new  arrangement. 
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INTRODUCTION 


This  memorandum  describes  the  program  6END  used  for  6 
component  force  and  moment  measurements  in  the  Transonic  Wind  Tunnel, 
and  it  is  intended  primarily  for  operators  of  this  wind  tunnel.  The 
program  is  for  use  with  strain  gauge  balances  using  first  order 
balance  interactions.  When  written,  the  8/1  installation  consisted 
of : 


DEC  PDP  8/1  with  extended  arithmatic  and  12K  of  core.  Plotter 
line  printer  LP3000,  Tektronix  Type  611  storage  display  and 
Tektronix  CIO  Polaroid  Camera,  Calcomp  Plotter  Type  565,  2 
minidiscs  type  DP32,  2  DEC  tape  transports,  multiplexer,  and 
ASR33  teletype.  This  has  some  time  ago,  changed  to  include 
24K  of  core,  ASR43  teletype,  high  speed  paper  tape  reader, 
dual  floppy  discs,  and  one  RK05  disc.  However,  the  original 
machine  dates  back  to  1968,  and  is  nearing  the  end  of  its  life, 
and  installation  is  scheduled  for  replacement  in  the  near 
future.  Nevertheless,  it  could  be  several  years  before  change¬ 
over  occurs,  and  meantime  the  present  operating  program  will  be 
used.  Further,  some  of  the  thought  behind  the  program  to  be 
discussed  here  may  be  useful  in  future  programs . 

The  program  6END  has  evolved  over  many  years  of  tunnel 
operation,  and  has  been  continually  modified  in  accordance  with 
tunnel  staff  wishes.  Operating  experience  showed  that  it  is  faster 
and  more  convenient  to  function  using  only  the  display  without  the 
plotter,  so  this  program  does  not  plot  on-line;  instead,  off-line 
plotting  and  cross-plotting  programs  were  developed,  and  are  now 
separate  programs.  Model  weight  corrections  have  been  incorporated, 
and  the  program  changed  to  suit  the  output  from  the  strain  gauge 
equipment.  The  program  has  been  in  use  for  many  years  and  has  been 
modified  from  time  to  time  in  response  to  users  recommendations. 

Listings  are  not  included,  since  the  program  occupies  3  fields  and 
is  too  big  to  include  here. 

2 .  PROGRAM  OPERATION 

Tables  must  be  compiled  as  in  Fig.  1  (see  Refs.  4,  5  and 
section  4) ,  keeping  in  mind  that  the  final  print  out  limits  coefficients 
including  C^a ,  to  a  maximum  of  99.  Tables  should  be  stored  on  disc 
as  on  ASCIifile  TABS,  and  all  signs  and  digits  should  be  entered, 
even  if  zero,  and  the  decimal  point  is  shown  by  the  dotted  line  in 
Fig.  1  and  should  not  be  entered.  Tabs,  carriage  returns  and  text 
may  be  entered,  but  will  not  be  saved  for  future  print  out.  It  is 
good  practice  to  store  the  disc  file  on  DEC  tape  so  that  it  can  be 
restored  to  the  disc  if  lost  from  there. 


F 
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The  program  uses  the  top  disc  for  data  storage  and  will 
check  to  see  that  the  top  disc  is  clear,  and  will  print  out  files  on 
the  top  disc,  but  it  is  more  efficient  to  clear  the  two  discs  of  all 
file  except  POP,  PIP,  .SYM,  TABS,  DET  (or  EDIT)  and  the  undeletable 
.CD  and  LOAD.  Using  PIP,  transfer  6END,  6011,  6012,  6111,  6211, 
and  print  out  of  disc  should  be  : 


FB=0463 

NAME  TYPE 

BLK 

8G 

POP  .SYS 

(1) 

0004 

PIP  .SYS 

(0) 

0025 

LOAD. SYS 

(0) 

0003 

.CD. .SYS 

(0) 

0006 

.SYM. SYS 

(0) 

0017 

DET'  .SYS 

(0) 

0025 

6011. SYS 

(0) 

0010 

6END.SYS 

(0) 

0037 

6012. SYS 

(0) 

0013 

6111. SYS 

(1) 

0034 

6211. SYS 

(2) 

0040 

TABS. ASCII 

0003 

*OPT- 

The  only  variation  permitted  is  that  the  number  of  blocks  for  TABS 
may  vary,  and  affect  the  quoted  number  of  free  blocks. 


With  the  computer  running  in  MONITOR,  type  6END,  carriage 
return,  and  the  program  will  load  and  stop,  having  found  and  loaded 
from  disc  6011  and  6211  (At  the  end  of  the  ident  and  tables  bit,  it 
loads  6012  and  6111) .  There  are  long  pauses  while  these  loadings 
occur.  It  is  good  practice  to  type  6END,  and  when  it  stops,  turn 
on  all  auxiliaries,  set  JOB  and  IDENT,  check  DEC  tapes  etc.  Then 
start  the  program  at  *201.  The  program  loads  6011  and  6211  and  then 
checks  that  the  top  disc  has  been  cleared.  If  not,  DISC  2  NOT 
CLEAR,  followed  by  the  names  of  the  offending  files  are  typed  out 
on  the  teletype  and  the  program  returns  to  MONITOR  to  allow  the  top 
disc  to  be  cleared. 


If  interface  power  has  not  been  turned  on  "NO  POWER"  appears 
on  the  teletype.  Turn  on  the  power  and  the  program  proceeds  normally, 
typing  DECTAPE? 


-  1 


The  question  is  whether  the  operator  wants  to  store  the 
ident  block  and  tables  on  DEC  tape  or  not.  For  verification  of 
tables  or  a  repeat  run,  or  for  preliminary  running,  DEC  tape  may 
not  be  required. 

Type  Y  for  Ves,  N  for  No  (no  CR  required).  If  Yes,  FIRST 
BLOCK  No  =  0  0  <t  0. 

CORRECT?  will  appear  and  requires  Y  or  N  as  before.  If  No,  FIRST  BLOCK 
Mo  *>  appears,  and  requires  4  digits  followed  by  C.R.  If  a  mistake 
in  typing  is  made,  type  RUBOUT  and  go  around  again.  JOB  NO  *  NN 
and  IDENT  NO  =  NNN  will  now  be  typed  out  on  the  teletype,  followed 
by  DATE  which  requires  NN  for  day  NN  for  month  NT!  for  year,  and  the 
program  inserts  spaces  between.  OPERATOR  will  accept  up  to  8  letters 
followed  by  CR  DETAILS  next  appears  and  still  expects  about  72 
characters  per  line  and  2  lines  maximum  -  i.e.  unchanged  from  ASR33. 

A  C.R.  is  required  to  terminate  each  line,  and  2  C.R.'s  in  all  must 
be  typed  for  the  program  to  proceed,  even  if  the  record  line  is  not 
used  for  text. 

It  is  essential  that  details  of  the  test  be  included  here, 
when  tunnel  runninq  begins,  and  that  it  be  stored  on  DEC  tape.  This 
is  the  only  means,  at  a  later  date,  of  identifying  the  data  that 
follows  on  the  DEC  tape.  Thus,  on  the  first  pass  of  the  program 
when  proper  data  are  expected  to  be  recorded,  sufficient  details 
must  be  included  to  identify  the  test  precisely.  Configuration 
details  should  be  included  and  enough  information  provided  so  that 
later  other  users  can  accurately  identify  the  data  -  otherwise,  the 
DEC  tape  shows  a  lot  of  data  that,  in  time,  may  be  meaningless. 

Later  passes  of  the  program  should  use  the  "Details"  section  to 
record  more  exact  details  of  the  data  to  follow.  In  all  tests,  the 
improtance  of  the  “Details"  cannot  be  overemphasized . 

By  now,  the  line  printer  has  printed  out 

JOB  NO  m  NN  IDENT  NO  =  NNN  Dates  NN  NN  NN 

OPERATOR:  ABC 

DETAILS:  The  Text  typed  in. 

The  tables  are  now  looked  for  as  "TABS"  on  the  discs,  and 
if  this  file  is  not  found,  the  program  expects  the  tables  from  paper 
tape  and  will  wait.  If  TABS  is  found,  or  paper  tape  fed  in,  the 
tables  are  read  into  Field  1,  printed  on  the  printer,  converted  to 
binary  and  stored  on  Field  2  starting  at  *201.  The  identification 
block  and  the  tables  are  now  stored  on  DEC  tape,  using  4  DEC  tape 
blocks  for  the  ident  segment  and  6  blocks  for  the  tables,  a  total  of 
10(10)  blocks,  provided  DEC  tape  has  been  requested  previously. 


If  DEC  tape  2  is  not  set  to  REMOTE,  or  WRITE  ENABLED,  or 
a  fault  occurs,  "DECTAPE  ERROR'-  will  appear  on  the  teletype,  and 
there  is  no  recourse  but  to  start  again.  The  program  now  loads 
6012  from  the  disc  into  core  as  an  overlay  into  Field  0.  Field  1 
previously  used  in  storing  tables  as  read  from  disc,  is  now  loaded 
with  6111  and  the  teletype  bell  rings  to  indicate  that  this  has 
been  achieved.  The  program  now  goes  to  *1001  or  Field  0. 

As  before  DECTAPE?  appears  and  if  Y,  FIRST  BLOCK  No  *  NNNN 

CORRECT? 

and  so  on  exactly  as  before.  However,  the  block  number  should  be 
correct,  unless  the  whole  program  has  been  reloaded.  The  DEC  tape 
query  here  simply  allows  the  operator  the  choice  of  doing  a  pre¬ 
liminary  run  with  full  print  out  but  no  DEC  tape  storage  and  is 
useful  for  checking  out  tables  etc. 

The  program  now  checks  to  see  if  the  display  is  turned  on, 
and  if  so  goes  to  Field  I.  If  desired,  typing  w  will  cause  the 
display  grid  to  be  drawn,  and  scales  adjusted  and  moved  (see  section 
3) .  The  program  now  waits,  apart  from  its  normal  computing  and 
displaying  Mach  number  on  the  control  desk,  until  the  "Record"  push 
button  is  pressed,  and  the  first  reading  must  be  initial  zero. 

While  more  than  one  set  of  initial  zeros  may  be  taken,  they  are  all 
printed  out  and  stored  on  DEC  tape,  which  tends  to  confuse  other 
programs  used  to  operate  on  the  data  later  on.  The  present  program 
singly  overwrites  initial  zeros  and  ignores  earlier  ones,  and  is 
unaffected. 

when  the  "Record"  P.B.  is  pressed  to  take  a  non-zero 
set  of  data  the  display  will  automatically  write.  All  display 
operations  are  now  permitted  until  the  final  zeros  are  taken.  Only 
one  set  of  final  zeros  is  possible  as  the  program  goes  off  to  compute 
all  the  data,  using  mean  zeros,  and  prints  it  out,  stores  it  on  DEC 
tape,  stores  it  on  the  discs,  and  cycles  back  ready  for  the  next  set 
of  data.  Test  block  number  is  printed  out  to  help  locate  the  data. 

At  present,  no  use  is  made  of  the  final  data  stored  on  the  discs. 
However,  with  increased  storage  capacity,  this  facility  could  be 
used  to  retain  data  for  several  configurations,  and  then  make  quick 
comparisons  using  the  display.  Such  programming  must  wait  until 
increased  storage  is  available. 

It  should  be  pointed  out  that  it  is  legitimate  to  restart 
the  program  at  *1001,  Field  0,  unless  the  program  has  been  destroyed 
by  some  means.  This  allows  the  tables  and  ident  part  of  the  program 
to  be  ignored,  and  the  tunnel  data  part  of  the  program  to  be  run. 
Therefore,  if  for  some  reason  a  tunnel  run  is  regarded  as  suspect, 
the  DEC  tape  block  number  can  be  set  to  the  right  value,  and  a  new 
run  carried  out. 
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DISPLAY  OPERATION 


3 . 1  Mode 


Normal  mode  of  operation  is  "Record’  mode,  and  this  is  the 
mode  automatically  entered  by  the  program.  Raw  data  are  printed  out 
and  stored  on  DEC  tape. 

"Look'-  mode  is  nrovided,  and  here  no  raw  data  are  printed 
or  written  on  DEC  tape.  This  mode  is  provided  so  that  the  tunnel 
operator  can  determine  the  shape  of  the  curves  and  decide  how  many 
and  where  data  points  are  needed.  For  this  purpose,  the  tunnel  may 
simply  be  accelerated  through  the  required  Mach  number  range  with 
the  operator  pressing  the  "Record"  first  button  as  needed  without 
bothering  to  obtain  precise  Mach  numbers .  Similarly,  at  constant 
Mach  number,  the  model  angle  of  attack  may  be  continuously  varied 
through  the  required  range  and  the  shape  of  the  curves  defined  by 
the  points  taken. 

Type  CTRL  R  for  Record  Mode. 

CTRL  L  for  Look  Mode. 

It  is  permissible  to  change  from  one  to  the  other  as  needed 
and  will  be  expected  by  the  program  any  time  after  initial  zeros 
and  before  final  zeros. 

3.2  Select  Abscissa 

Type  V  A  then  W  to  display  coefficients  versus  a. 

Type  V  B  then  w  to  display  coefficients  versus  g. 

Type  V  M  then  W  to  display  coefficients  \ersus  Mach  number. 

The  program  loads  in  the  V  M  state.  It  is  permissible  to 

go  from  one  to  another  and  back  again. 

3.3  Write  Instruction 

Type  W.  This  will  renew  the  display,  or  write  the  new  one 
requested  in  (2)  above. 

3 . 4  Deletions 

Type  L  N 
M 

C  . 

/  followed  by  w. 


3 . 6  Movim 


Type  t  to  move  up. 

Type  !  to  move  down. 

Type  >  to  move  right. 

Type  <  to  move  left. 

Type  RUBOUT  to  centre  both. 

Each  time  one  of  these  is  typed,  the  display  is  moved  by  one  grid 
line  on  the  display,  whatever  that  is  at  the  time  of  moving.  It 
is  thus  faster  to  move  the  display  first  and  then  expand,  then  the 
other  way  round. 

Initial  scales  are  : 

M:#  0.4  to  1.4  in  increments  per  grid  line  of  0.2. 

a:#  0  to  60  deg.  in  increments  per  grid  line  of  10  deg. 

3:#  0  to  60  deg.  in  increments  per  grid  line  of  10  deg. 

Coefficients:-  -4  to  +4  in  increments  per  grid  line  of  1. 

During  moving  or  zooming  points  may  be  lost  from  the 
display.  Such  points  are  never  truly  lost  and  will  reappear  when 
the  display  is  suitably  moved  or  scaled.  (In  part,  the  display  has 
a  visible  window  in  a  larger  storage  area  and  program  protection 
takes  over  if  the  larger  area  is  exceeded) . 

It  is  also  quite  possible  to  get  a  mixture  of  zooming  and 
moving  such  that  it  becomes  most  confusing  to  operate  further  and 
in  such  cases  it  is  frequently  quicker  to  type  vmw  (or  similar)  and 
start  over  from  the  beginning. 

3.7  XI  or  X10  Facility 

When  the  display  is  operated,  it  will  be  seen  that  some 
coefficients  are  plotted  X10.  Any  coefficient  may  be  displayed 
XI  or  X10  by  changing  the  contexts  of  locations  #3000  to  #3005  on 
Field  1.  0  in  these  locations  produces  Xlu,  7777  produces  XI. 

Locations  are  : 

#3000  L  3003  C 

3001  K  3004  N 

3002  D  3005  1 

To  do  this,  after  the  Field  1  program  has  been  loaded  and 
before  taking  initial  zeros,  stop  the  PDP8/I.  Toggle  the  required 
values  in  to  the  appropriate  locations  and  restart  at  #1001  or  Field  0. 


This  program  change  will  be  lost  if  the  program  is  run  from  the 
beginning  again,  and  to  avoid  frequent  toggling  it  is  better  to 
SAVE  the  toggled  urogram  as  6111. 


3.8  Restrictions 

Pressing  the  'Record11  P.B.  causes  data  to  be  written  on 
the  discs,  operating  the  display  requires  data  to  be  read  from  the 
discs.  Efforts  to  provide  protection  have  not  been  successful, 
and  it  is  therefore  best  not  to  push  the  P.B.  while  the  display  is 
writing  or  use  the  teletype  to  issue  display  instructions  while  the 
recording  program,  which  includes  redrawing  the  display,  is  proceeding 
In  practice,  this  is  no  hardship  as  the  operator  normally  presses 
the  button  to  get  his  reading  and  then  looks  at  the  display,  makes 
any  changes  needed  to  this,  and  then  takes  the  next  reading. 

As  at  present  programmed,  as  soon  as  final  zeros  are 
recorded,  the  display  becomes  inoperative.  Although  the  actual 
picture  remains,  the  program  will  not  listen  to  the  teletype  because 
the  program  has  gone  on  and  is  ready  for  the  next  run.  Consequently, 
if  the  display  is  to  be  photographed,  the  display  must  be  set  up 
to  the  form  needed  before  taking  final  zeros. 

3 . 9  Photographing  the  Display 

A  Polaroid  camera  is  available  to  take  photographs  of 
the  display,  and  gives  a  print  1/2  full  size.  Remove  both  lens  caps, 
press  camera  against  screen,  and  press  shutter  and  follow  the 
Polaroid  directions  is  all  that  is  normally  required.  However,  the 
screen  blanks  off  for  a  short  time  every  90  seconds  approximately, 
and  so,  if  possible  type  W  just  before  taking  the  photo,  or  if  this 
is  not  possible,  wait  till  the  display  comes  on  after  blanking  out, 
and  then  photograph. 

3.10  Malfunction 


Every  so  often,  about  once  in  2  or  3  days,  the  display 
will  draw  only  part  of  the  picture,  and  the  program  locks  up.  Press 
"ERASE"  and  the  program  v/ill  continue.  The  fault  appears  to  be  a 
hardware  fault  in  the  interrupt  circuitry,  but  is  too  intermittent 
to  make  fault  finding  practical. 

4 .  DESCRIPTION  OF  PROGRAM 

4 . 1  General 

This  program,  6END,  occupies  all  fields  and  field  0  uses 
one  overlay.  Terminology  used  for  the  component  parts  of  the  program 
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6  for  6  six  component  program 
for  field 

I  for  insert 

N?  for  number  of  the  insert. 

Thus  field  2  is  6211,  field  1  is  6111,  and  field  0  uses  6011  and 
6012.  After  6END  is  loaded,  the  program  is  self  loading,  and  it 
is  therefore  easy  to  add  future  inserts  or  overlays  as  needed.  Thus 
6112  could  be  a  plot  program  to  overlay  the  display  program  and 
provide  plotted  data  at  the  end  of  the  tunnel  run. 

6 END,  the  first  part  of  the  program  loaded,  contains  the 
subroutine  section-*  which  is  common  to  most  tunnel  running  programs, 
and  is  therefore  kept  as  a  separate  entity.  It  contains  B.C.D.  to 
binary,  binary  to  B.C.D. ,  basic  multiply  and  divide,  Reynolds  number 
and  Mach  number  subroutines.  These  are  largely  the  same  as  in 
Ref.  3.  All  are  obvious  except  the  Mach  number  subroutine  and  its 
operation.  Wherever  6END  is  running  and  the  interrupt  is  turned  ON, 
it  is  interrupted  about  twice  a  second  by  pulses  from  an  electronic 
clock.  When  this  interrupt  occurs,  tunnel  static  and  total  pressure 
which  are  continuously  measured  by  automatic  manometers,  are  read 
in,  Mach  number  computed  and  then  displayed  on  the  tunnel  control 
desk.  Thus,  all  the  time  6EMD  is  performing  other  functions,  it  is 
still  producing  tunnel  Mach  number  for  the  plant  operator. 

In  this  connection  it  should  be  noted  that  6 END,  which  uses 
the  same  Mach  number  subroutine,  computes  Ma  and  uses  this  to  compute 
coefficients.  The  Mach  number  printed  out  by  the  program  should  not 
be  squared  -  it  is  much  more  accurate  to  get  M3  from  static  and  total 
pressure  -  for  any  further  computations. 

These  subroutines  occupy,  on  Field  0,  *50-*75,  *106~*131, 
*6244-*7551,  the  page  0  locations  being  used  for  addressing  the 
subroutines . 

The  Identification  section  is  much  the  same  as  before, 
but  tables  must  now  be  in  the  form  shown  in  Fig.  1,  and  are  expected 
to  be  on  disc  as  on  ASCII  file  TABS.  To  eliminate  the  weight  correction 
terms,  zeros  should  be  entered,  and  to  operate  with  less  than  6 
components,  the  strain  gauge  balance  tables  should  contain  zeros  in 
the  appropriate  places.  To  produce  the  inverse  balance  parameters 
matrix  see  Refs.  4  and  5.  The  Identification  section  is  actually  now 

6011,  occupying  *1000-*2777. 

The  record  part  of  the  program  occupies  *1000-*3777  and  is 

6012.  Thi3  part  of  the  program  sorts  out  zeros  from  readings,  types 


headings,  checks  whether  the  display  is  on  or  not,  sends  the  program 
off  to  compute  etc.  When  the  "Record'  push  button  is  pressed,  various 
delays  are  initialized,  all  6  A  to  D  converters  in  the  strain  gauge 
equipment  are  started,  and  the  digital  thermometer  solenoid  energized. 
After  an  appropriate  time  delay  an  interrupt  is  issued,  and  the 
program  identifies  the  interrupt  and  goes  to  the  actual  "Record'1 
subroutine,  which  starts  at  *1055.  The  push  button  is  inhibited  for 
a  period  of  about  5  seconds  to  prevent  multiple  readings  if  it  is 
dropped,  and  to  ensure  that  only  1  set  of  data  is  taken  at  a  time. 

The  display  program  resides  in  Field  1,  and  is  6111, 
occupying  *1000-*7777.  Field  2  contains  the  main  computing  section, 
and  stores  zeros,  tables  etc  It  is  6211,  occupying  the  whole 
field. 

Computation  is  based  on  Ref.  5,  so  that  using  the  notation 
and  definition  given  in  Appendix  1. 

Z  ■  k1„„  S„  (r  -  r„  )  +  -  +  K1  ST  (rr  -  r.  ) 

ZZ  Z  Z  Zo  ZL  L  L  Lo 

with  similar  expressions  for  M,X,Y,M,L. 

n  =  xz.z  +  x^.m. 

♦  -  xy.y  +  an.n 

♦*-  XL.b 


sin  a  =  sin0  cosn  cos<j>  +  cos0  sing 
sin  8  »  sin8  cosi//  sin$  -  cos0  sinty 
sin  a1  =  sina  secB 


2.0  36  Z 

CZ  ®  0.700  p  M  3  S 

o 

with  similar  expressions  for  C  ,  C 

A  Z 


„  _  2 .036  M 

"M  =  0.700  p  f  J  Sc 

o 

-  .  <PB  '  P>  SB 

"X„  0.700  p  f:  a  S 


-IS¬ 


C'  =  c  -  c 

X  X  X 


B 


2.036  N 


n  0.7000  p  M  a  Sb 
o 

and  a  similar  relation  for  C„ 


C  =  -C  secB  /cos3  a  -  sin3  B  +  C'  sina  secB 
Li  Z  X 


C  =  -C '  /cos3  a  -  sin3  B  -  C  sinB  -  C  sina 
U  X  Y  z 


C  =  C '  tanB  /cos3  a  -  sin3  B  -  C  cosB  +  C  sina  tanB 

C  X  Y  Z 


The  arcsine  subroutine  used  is  based  on 


arcsin  X  =  ^  -  /l  -  X  ¥(X) 


where 

Y(X> 

as 

a0  +  ax  X  +  a2 

where 

a 

o 

= 

1.5707,  8786 

al 

= 

-  .2141,  2453 

a2 

= 

+  .0846,  6649 

a3 

= 

-  .0357,  5663 

a4 

= 

.0086,  4884. 

This  subroutine  gives  X  to  suit  the  size  subroutine  used,  with  X  in 
degrees  X100,  and  shifted  left  9. 

The  double  precision  square  root  subroutine  is  based  on  ; 

xi«  -  1/2 

It  should  be  noted  that  the  relation  used,  and  the  arcsin  subroutine, 
limit  the  program  to  -90°  <  0  <  +90°. 

4 . 2  Weight  Corrections 

A  conventional  strain  gauge  balance  measures  the  total  load 
on  the  model,  and  the  total  load  consists  of  the  wind  load  and  the 
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weight  of  the  model.  The  original  program  ignored  the  weight  of 
the  model,  which  meant  that  operation  at  fixed  model  attitude  and 
varying  Mach  number  was  necessary,  and  model  weight  is  then  eliminated. 

The  program  modifications  are  based  on  equations  derived 
by  D.A.  Secomb. 

The  model  attitude  relative  to  horizontal  is  computed  from  : 

6  =  9  +  n  eos<t>  -  V  sin<J>  -  6 

3  O 

0  =  0  +  0'  -  0 

a  o 

where  6  and  0  are  the  nominal  pitch  and  roll  angles 
recorded  for  each  data  point. 

6  and  0  are  the  values  of  the  nominal  pitch  and  roll 
angles  for  which  tiie  model  reference  system  is  horizontal,  often 
zero. 


n,  ¥  and  0'  are  the  small  angular  deflections  due  to  the 
combined  effects  of  the  aerodynamic  load  and  the  change  in  model 
attitude . 


The  weight  components  are  given  by  : 

Z  =  -W  cos  0  cos  0 
W  a  a 

M_.  =  -W  (Z  sin  0  +  X  cos  0  cos  0  ) 

tv  g  a  g  a  a 

X  =  -W  sin  0 

W  a 

Y  =  W  cos  0  sin  0 

W  a  a 

ML,  =  W  (X  cos  0  sin  0  +  Y  sin  0  ) 

W  g  a  a  g  a 

L  =  W  (Y  cos  0  cos  0  -  Z  cos  0  sin  0  ) . 

W  g  a  a  g  a  a 

where  W  is  the  effective  weight  of  the  model,  and  (X  ,  Y  ,  Z  )  are 
the  coordinates  of  the  centre  of  gravity  relative  to9thegmod£l 
reference  system,  which  is  the  axis  system  to  which  the  balance  matrix 
and  stiffness  factors  are  referred. 

These  equations  are  computed  in  the  actual  program  as  : 

Zr.  =  (V?  cos  6  )  cos  0 

w  ct  d 

\  =  \Zg-ZWXg 

X  *  -W  sin  0 

w  a 

Y  =  (W  cos  9  )  sin  0 

”  a  a 


g 


These  values  are  computed  for  wind-off  zeros,  and  stored,  and  then 
computed  for  each  data  point,  and  the  increments  (etc.)  subtracted 
from  the  measured  2  (etc.)  to  give  the  aerodynamic  component . 

The  validity  of  the  tables  may  be  checked  as  follows  : 

(i)  Take  initial  wind  off  zeros  at  a  given  model  altitude. 

(ii)  Set  p  and  H  to  give  a  fake  Mach  number,  and  this  should 
remain  constant. 

(iii)  Change  the  model  altitude  setting  and  take  data  points  as 
required,  still  with  wind  off. 

(iv)  Return  Mach  number  to  zero,  model  altitude  to  its  initial 
value  and  take  final  zeros. 

The  program  will  now  compute  coefficients  and  these  should 
all  be  zero  within  the  limits  of  resolution  of  the  whole 
system.  If  not,  one  or  more  of  the  values  entered  in  the 
tables  is  incorrect  and  needs  modifying. 

4 . 3  Storage  Information 


4.3.1  Core 


1 

Name  j 

Starting 

Location 

Tables 

2 

201 

Raw  data  First  zero 

2 

551 

Final  zero 

2 

1001 

Reading 

2 

1301 

i 

j  Computed  data 

2 

1301 

j  Mean  zero 

2 

1400 

Wind  off  weight  correction  terms 

i 

2 

1505 

Only  one  data  block  is  used  and  computed  data  is  written 
over  the  raw  data. 


Zeros 

Disc  Starting 

Blocks 

Used 

Address 

Available  (dec . ) 

Look 

Computed 

Initial 

0010,1003 

410 

Record 

Raw 

- 

0014,0200 

208 

Computed 

Initial 

0016,0100 

2oe 

Computed 

Mean 

0016,0100 

208 

|  Top  of  Disc  2 

) 

is  0017,7775 

No  protection  is  provided  in  the  program,  and  no  indication 
that  the  number  of  blocks  provided  has  been  exceeded.  Therefore,  it 
is  quite  permissible ,  to  use  the  whole  826  blocks  in  Look  mode,  but 
these  data  will  be  written  over  when  Record  mode  is  entered.  If  a 
very  large  number  of  Record  blocks  is  required,  care  should  be  taken 
that  the  raw  data  do  not  write  over  the  computed  data  -  which  will 
be  immediately  obvious  on  the  display.  Otherwise,  appropriate 
changes  to  the  DISSET  subroutine,  *2205  Field  0,  may  be  made  to 
double  the  Record  storage  and  delete  the  Look  storage.  This  involves 
changing  the  raw  data  address  to  the  Look  disc  address,  and  the 
computed  data  address  to  the  original  raw  data  address,  which  is 
easier  to  do  than  to  describe. 

In  practice,  a  typical  tunnel  run  uses  around  20  -  50 
blocks  and  so  very  ample  storage  is  already  provided. 

4,3.3  DEC  tape 

According  to  DEC,  DEC  tape  storage  is  given  by  : 

N  =  212112  -  2. 

(N  +  15) 


where  N  is  number  of  words  per  block  and  must  be  divisable  by  3. 

N  is  number  of  blocks  per  reel  of  DEC  tape,  standard 
block  size  in  this  program  is  39  words,  giving  3926  blocks  per 
reel  of  DEC  tape. 


As  a  general  principle,  if  raw  data  are  to  be  stored,  it 
should  be  stored  as  soon  as  possible  and  as  acquired.  Therefore, 

6END  writes  the  raw  data  in  the  format  shown  on  Fig.  5  on  DEC  tape 
as  soon  as  possible,  and  then  proceeds  to  print  it  out  and  conduct 
further  manipulations. 

Raw  data  format  for  the  disc  storage  is  shown  in  Fig.  6 
and  is  now  all  binary,  except  for  SERIAL  number,  and  the  numbers 
not  over  written  as  indicated  in  Fig.  6. 

The  format  for  computed  data  on  both  disc  and  DEC  tape, 
is  that  shown  in  Fig.  5. 

4 . 4  Protection  and  Limitations 

1.  Strain  gauge  balance  outputs  cannot  exceed  9999,  but  in 
practice,  the  balance  must  not  be  overloaded.  Set  the  2 
high  order  digits  on  the  thumbwheel  switches  on  the  strain 
gauge  equipment  for  all  channels  to  be  used,  to  a  value 
corresponding  to  maximum  permitted  outputs  for  the  balance 
being  used.  A  beeper  will  sound  and  the  LED  light  up  on 
the  offending  channel  if  the  preset  value  is  reached. 

2.  Static,  total  and  base  pressure  are  restricted  to  one 
word  -  about  42.00. 

3.  Reference  length  $  16  inches. 

4.  M  <  0.10  is  treated  as  a  zero,  but  in  normal  operation, 

M  =  0.000  is  required. 

5.  For  8  >  90°,  the  computing  program  will  give  incorrect 
answers;  raw  data  is  unaffected,  but  the  maximum  set 
table  value  is  presently  99.99°. 

6.  Poll  angle  can  be  up  to  °  200.0°. 

7.  DEC  tape  must  be  used  for  data;  format  must  be  39(10) 
words  per  block. 

8.  For  Look  mode,  410  points  are  allowed;  in  Record  mode, 

208  points  (see  4.1). 

9.  Do  not  operate  Record  push  button  while  display  is  writing. 

10.  Do  not  issue  teletype  instructions  while  display  is  writing 
or  recording  of  data  is  proceeding. 


5. 


conclusion 


The  current  version  of  the  program  used  for  six  component 
measurements  using  strain  gauge  balances  and  the  DSC  8/1  installation 
in  the  Transonic  Tunnel  has  been  described.  Details  of  the  display 
programs  and  other  changes  to  early  programs  are  also  described, 
as  well  as  operating  details.  The  program  described  has  been  in 
use  for  many  years,  being  modified  from  time  to  time  to  suit  users 
recommendations . 
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APPENDIX  I 


r 
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NOTATION  AMP  DEFINITION 
(See  Reference  5) 


I 

1 


1.1 


7 

M 

X 

Y 

N 

L 


Body  Axes  Force  System 

Normal  Force,  positive  downwards. 

Pitching  Moment,  positive  nose-up. 

Axial  Force,  positive  up-wind. 

Side  Force,  positive  to  starboard. 

Yawing  Moment ,  positive  nose  to  starboard 

Rolling  Moment,  positive  for  drop  of  starboard 
wing. 


±  xxx. xxx  lb. 

±  xxx. xxx  lb. in. 
±  xxx. xxx  lb. 

±  xxx. xxx  lb. 

±  xxx. xxx  lb. in. 
i  xxx. xxx  lb. in. 


I 


1 


1 . 2  Attitude  Angles 

6  Sting  root  or  nominal  pitch  angle,  positive  ±  xx.xx  deg 

nose-up . 

Sting  root  or  nominal  roll  angle,  positive  for.  i  xxx.t  deg. 
drop  of  starboard  wing. 

4>m  Model  roll  angle,  wind-off,  positive  for  drop  t  xxx.x  deg. 

of  starboard  wing. 

a  Angle  of  incidence  :  i.e.  angle  between  the  i  xx.xx  deg. 

wind  vector  and  its  projection  on  the  chord 
plane,  positive  nose-up*. 


0  Angle  of  sideslip  :  i.e.  angle  between  the  wind  ±  xx.xx  deg. 

vector  and  its  projection  on  the  model  plane  of 
symmetry,  positive  nose  to  port*. 


Alternative  angle  of  incidence  :  i.e.  angle  ±  xx.xx  deg. 

between  the  model  axis  and  the  projection  of 
the  wind  vector  on  the  model  plane  of 
symmetry* . 


* 


* 


These  definitions  apply  strictly  only  in  the  range  -90°  to  +90°. 


Anqular  deflection  under  load  of  model,  ±  x.xx  deg. 

neasured  in  the  plane  of  symmetry,  positive 

nose-up. 

Angular  deflection  under  load  of  model ,  ±  x.xx  deg. 

measured  perpendicular  to  the  plane  of  symmetry, 
positive  nose  to  starboard. 

Angular  deflection  under  load  of  model  measured  ±  x.xx  deg. 

around  the  axis  of  symmetry,  positive  for  drop 
of  starboard  wind. 

Model  roll  angle,  wind  on,  positive  for  drop  ±  xxx.x  deg. 
of  starboard  wing . 

(♦-♦'  +  <J>B) 

Sting  Balance  Stiffness  Parameters 

deflection  per  degree  by  Z  load  ±  .xxxx  deg. /lb. 

deflection  per  degree  by  M  load  ±  .xxxx  deg. /lb. 

in. 

deflection  per  degree  by  Y  load  ±  .xxxx  deg. /lb. 

deflection  per  decree  by  W  load  t  .xxxx  deg. /lb. 

in. 

deflection  per  degree  by  L  load  ±  .xxxx  deg. /lb. 

in. 


Directions  of  Hind  Axes 

The  lift  direction  is  that  perpendicular  to  the 
wind  vector,  in  the  plane  of  symmetry,  positive 
upwards  for  a  model  in  its  datur  orientation. 

The  drag  direction  is  parallel  with  the  wind 
vector,  positive  downwind. 

The  cross-wind  direction  is  perpendicular  to 
each  of  the  above  directions,  and  is  positive 
to  port  for  a  model  in  its  datum  orientation. 


1 


c 

; 

1 


< 


i 


I 


i 


AI.3 


.5 

Wind  Axes  Forces 

L 

Lift  force 

± 

XXX. XXX 

lb. 

D 

Drag  force 

± 

XXX. XXX 

lb. 

C 

Cross-»/in<l  force 

± 

XXX. XXX 

lb. 

are  the  resolutes  of  the  total  aerodynamic 
force  in  the  defined  directions. 


L .  6 

Model  Geometry 

S 

'•'ing  area 

XXX. XXX 

2 

xn. 

SB 

Base  area 

XX .XXX 

2 

xn. 

c 

Reference  chord  length 

XX .XXX 

in. 

b 

Reference  span  lenqth 

XX. XXX 

in. 

*R 

Reference  length  for  Reynolds  Number 

XX. XXX 

in. 

L.7 

Aerodynamic  Quantities 

P 

Reference  static  pressure. 

XX. XX 

in . Hg . 

H 

Reference  total  pressure . 

XX. XX 

in.Hg. 

T 

o 

Total  temperature. 

XX. X 

deg.C. 

PB 

Model  base  pressure. 

XX. XX 

in.Hg. 

M 

o 

Mach  number  based  on  p  and  H. 

x.xxx 

0.7 

pM  ^  Kinetic  pressure, 
o 

XX  .  XX 

in.Hg. 

Re 

Reynolds  number. 

XX. XXX 

million 

2 

(The  conversion  factor  from  in. Ha.  to  lb. /in.  is 
2.036) 


x.xxxx 


AI.4 


1 . 8  Body  Axes  Force  Coefficients 

Normal  Force  Coefficient. 

C  Pitchino  foment  Coefficient, 

m 

Cx  Axial  Force  Coefficient. 

Cy  Side  Force  Coefficient. 

C  Yawing  i'oment  Coefficient, 

n 

C^  Rolling  Moment  Coefficient. 

C  Base  Force  Coefficient. 

B 

Cx*  Corrected  axial  force  coefficient. 


±  X . xxxx 
±  X  .xxxx 
±  X  .  xxxx 
i  x.xxxx 
t  x.xxxx 
±  x.xxxx 
±  x.xxxx 


i • 9  rind  Axes  Force  Coefficients 

CL  Lift  Coefficient. 

CQ  Drag  Coefficient. 

Cc  Cross  Wing  Force  Coefficient 

For  all  coefficients  greater  than  9.9999  print 
out  is  t  xx. xxx 


±  X . xxxx 
±  x.xxxx 
±  X . xxxx 


1.1C  Bridge  Scale Readings 
rz  Reading  of  7  scale, 
r^  Reading  of  M  scale 
rx  Reading  of  X  scale 
ry  Reading  of  Y  scale. 
rN  Reading  of  M  scale. 


i  xxxx 
t  xxxx 
±  xxxx 
i  xxxx 
±  xxxx 


r, 


Reading  of  L  scale. 


±  xxxx 


Suffix 


o 


indicates  wine  off  zero  e.g.  r 
o  Z 

o 

Bar  indicates  mean  of  initial  and  final  values 
e-g-  r„ 


1.11  Bridge  Sensitivity  Factors 
Sp  Sensitivity  of  P  channel 

where  P  may  be  Z,  M,  X,  Y,  M,  L. 

Load  for  full  scale  deflection  of  9999. 


±  XXX. XXX 


1.12  Inverse  Balance  Parameters 


K',  K*  etc.  (see  text). 
ZZ  ZM 


1 , 13  Weight  Correction  Terns 
0 

o 
*0 

w 


Pitch  angle  when  model  is  horizontal 

± 

XX. XX 

deg 

Roll  angle  when  model  is  horizontal 

t 

XXX. X 

deg 

Effective  model  weight 

± 

XXX. XXX 

lb. 

Axial  coordinate  of  C.G. 

± 

XX. XXX 

in. 

Lateral  coordinate  of  C.G. 

± 

x.xxx 

in. 

Vertical  coordinate  of  C.G. 

± 

x.xxx 

in. 
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16  Abstract 

Thin  memorandum  describes  the  current  version  of  the 
program  used  for  6  component  measurements  using  strain  gauge 
balances  and  the  Pi, PE/I  instillation  in  the  Transonic  Wind  Tunnel. 
The  program  is  an  updated  version  of  ARL  Tech.  Memo.  A2641  and  has 
been  in  u3e  fer  several  years.  Changes  to  the  original  program 
include  on-line  display,  model  weight  correction,  and  operation  with 
the  Pollock  strain  gauge  equipment2.  Although  the  8/1  is  scheduled 
for  replacement,  it  may  be  some  time  before  the  changeover  to  a  new 
installation  can  occur.  Also,  some  of  the  program  detail  may  be 
useful  to  the  new  arrangement. 
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